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Large cutoff effects of dynamical Wilson fermions * 
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We present and discuss results for cutoff effects in the PCAC masses and the mass dependence of ro for full 
QCD and various fermion actions. Our discussion of how one computes mass dependences - here of ro - is also 
relevant for comparisons with chiral perturbation theory. 



1. INTRODUCTION 

In the quenched approximation, large cutoff- 
effects were observed with the original Wil- 
son action. In particular, they appeared as 
a strong dependence of the PCAC-mass, m = 
^{P\df^A^\a) / {(3\P\a), on the external states, 
I a), 1/3), a dependence which has to be absent in 
the continuum limit pp. These a-effects could be 
reduced to a tolerable level by non-perturbative 
0(a)-improvement 0. The necessary improve- 
ment coefficient Csw is now known also for the Wil- 
son gauge action with A'^f = 2 flavors of dynamical 
Wilson fermions (action "W/SW") and with 
the Iwasaki gauge action with A^f = 3 (I/SW) 
It was pointed out already in [3] that O(a^) effects 
may be sizeable at the lowest value of (3 consid- 
ered, which roughly corresponds to a ~ 0.1 fm. 
Here we report on O(a^) effects with dynamical 
fermions, for W/SW and I/SW actions and for 
comparison we consider also the original Wilson 
(W/W) and Kogut Susskind (W/KS) actions. 



2. PCAC MASSES 

A family of external states, |a), \(3), can be re- 
alized in the Schrodinger Functional (SF) ^E]. 
The (bare) PCAC-mass, 

_ dafAixo) + CAad^dofpixo) 



2/p(xo) 



(1) 



xo=L/2 



then depends on the resolution L/aoi the x L 
space-time, on the spatial periodicity angle 9 of 
the fermion fields and the Dirichlet boundary con- 
ditions for the gauge fields. In eq. ^ the stan- 
dard SF correlation functions of the axial vector 
and the pseudoscalar density /a, /p enter. Chiral 
symmetry predicts that m is (apart from cutoff ef- 
fects) independent of L/a as well as of the other 
parameters, as long as the gauge coupling and the 
bare quark mass are kept constant. Denoting by 
TTioo the PCAC mass computed in large volume 
with \a) the p = one tt state and \/3) = |0), one 
expects in particular 
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in the 0(a) improved cases. For W/W this differ- 
ence is instead 0(a). Indeed, after 0(a) improve- 
ment Am is not too large for a « (2GcV)^^ « 
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Am[MeV] x ^^ij^ 
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12(1) 


w/sw 


5.2 
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28(1) 


w/w 


5.5 
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96(4) 


I/SW 


2.2 
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-2(4) 



Table 1 

Am at a^^ ~ 2GeV, L/a — 8. Values for moo have 
been taken from EIIZIIH]- Except for Nf = 0, the 1- 
loop expressions for the improvement coefficient ca 
are used. For the first three lines we have 9=1/2 
and no background field in the SF, while for the last 
line, 6 = and the background field of (2) was chosen. 



O.lfm in the quenched approximation. However 
for A^f = 2 it is more than a factor two larger 
(Table and, although not shown here, this 
roughly holds also for other such mass differ- 
ences. Considering several other differences, we 
convinced ourselves that the only perturbatively 
known value for ca is not the origin of the large 
cutoff effects seen for W/SW and iVf = 2. 

We then considered the Iwasaki gauge action 
together with the values of Csw used in (Hj- We 
found a much smaller Am, see Tabled 

3. MASS DEPENDENCE OF ro 

In lattice gauge theory computations, the de- 
pendence of ro on the mass(es) of the dynamical 
quarks has been neglected so far. However, in 
contrast to the qq force at very short distance r 
where only an effect of relative size ©(ofgrn^r^) is 
present, for the force at r w 0.5 fm this depen- 
dence is not obviously that small. It should be 
computed by lattice gauge theory. Due to spon- 
taneous chiral symmetry breaking, a linear term 
in m is expected for small quark masses. 

In order to compute the mass dependence on 
the lattice, one must not introduce spurious mass 
dependences through renormalization. A mass in- 
dependent renormalization scheme has to be cho- 
sen. In perturbation theory an example is pro- 
vided by the lattice minimal subtraction scheme 
with renormalized coupling 



= {1 ^ 2boHa^i)gl + . . .} gl . 

The above equation shows that the lattice spac- 
ing is fixed by go, independently of the quark 
mass. Apart from a small 0(a)-correction in 
0(a)-improved QCD, which is due to [S] 



5o = (1 + &gamq)go : 



(4) 



the lattice spacing is determined by go alone in 
a mass-independent renormalization scheme. We 
note in passing that it is also mandatory to choose 
such a scheme in a comparison of lattice QCD re- 
sults with chiral perturbation theory predictions. 

To define the mass-dependence, we first intro- 
duce a reference value mrcf of the quark mass via 
rpTTipglm^mret — which we later put to i^T = 3 
in our numerical evaluation. Other choices are 
possible. The physical mass dependence of ro 
may then be described by the function pcont (x) = 
ro(m)/ro(mrof) where x — m/m-cci is finite since 
in the corresponding renormalized ratio the com- 
mon Z-factor (defined in the massless scheme) 
cancels. A precise definition of Pcont is provided 
by (ro = ro/a) 

Pcont{x) ^ lim p(x, l/fo(mi.cf)) , (5) 

l/fo(mrcf)^0 



p{x, l/fo(mi.ct)) = 



fo(m) 



ro(mrcf) 



(6) 



go fixed 



As a first step we try to understand whether the 
mass dependence is a large effect and quantify 
scaling violations. Since we wanted to use results 
from the literature, it was easier to consider 
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9ut{p) = Zg{gQ,afi)gl 
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which agrees with p at small quark masses and 
lattice spacings due to to ~ "ipg. In practice 
only go could be kept constant in our evaluation 
of i?, assuming that the bganiq term is negligible 
(6g « 0.02 + O(g^) for W/SW 0). It turns out 
that R is well represented by the Taylor expansion 

R{X, l/fo(TOref)) = 1 + s(l/fo(TOref)) {x - 1) (7) 

for the data considered [H1E1IZIEI|> namely 
I/SW x G [0.4, 2.0] , fo(mrcf) e [1.8, 5.3] 
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Figure 1. Lattice spacing dependence of slopes s. 

W/W X G [0.8, 1.6] , fo(m,.cf) « 4.6 
W/SW X e [0.6, 2.0] , fo(TO,.cf) « 4.8 
W/KS X e [0.3, 2.0] , fo(TOi.cf) e [1.8, 4.9] 

Our figure shows the available results for the slope 
s. Both the I/SW numbers and the W/KS ones 
point to a very small |s| in the continuum limit 
with a rather strong a-dependence and unclear 
mutual consistency. The two other actions appear 
to show even larger cutoff effects at a « 0.1 fm. 

4. DISCUSSION 

Our results are not easy to interpret, apart 
from the general observation that dynamical 
fermions may introduce significant cutoff effects 
at a 0.1 fm. We recall some other findings. 

• UKQCD found a very small and strongly 
mass dependent 0+"*" glueball mass at a « 
0.1 fm with Nf = 2 W/SW [U]. 

• For the same action and parameters, the 
auto-correlation times increase as the quark 
mass is increased. 

• The pure SU(3) gauge theory has a phase 
transition in the (3-p, Pa > half-plane not 
very far from a ~ 0.1 fm. 

• The JLQCD collaboration found a phase 
transition with Nf = 3 W/SW pSI, pre- 
venting simulations unless a is quite small. 

These hints lead us to the conjecture that for 
N{ = 2, W/SW, a « 0.1 fm, one may be close 
to a phase transition and thus suffer from cut- 
off effects which are not necessarily described by 
the Symanzik expansion. Concerning the PCAC 



masses (see section 2), I/SW is much better. 
However, we are lacking an understanding why 
this is so and our figure shows significant a-cffects 
in the mass dependence of tq for I/SW. 

At the present time we draw two conclusions. 

1) The issue of what is a good action for full lat- 
tice QCD deserves more attention and one should 
investigate scaling properties of actions which are 
used to compute physical observables. 

2) The mass dependence of tq and in particular 
the lattice artifacts contained in this mass depen- 
dence has to be understood before drawing con- 
clusions on the chiral behavior of quantities such 
as and rn^/ra. Section 3 outlines what needs 
to be done but our numerical investigation should 
be considered only as a first step. 
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